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14  Abstract 


II.  OntnkitiM  Slati 


5Avcj  Corporation  *s  Systems  Division  designed  and  developed  an  engineered 
Pulsed  Acoustic  Vo"tex  Sensing  System  (PAVSS).  This  system  Is  capable  of 
real-time  detection,  tracking,  reco  -ding,  and  graphic  display  of  aircraft 
trailing  vortices. 


Tt is  volume  of  the  report  presents  the  results  of  two  subcontractor 
studies  directed  toward  development  of  improved  vortex  tracking  software 
techniq  es  for  the  PAVSS.  The  volume  recommends  the  incorporation  of 
several  improvements  in  the  software.  The  subcontractor  final  reports 
(Scope  Electronics,  Inc.  and  Arcon  Corporation)  are  furnished  as  appendixe, 
to  this  volume 

Other  volumes  in  this  final  report  are  as  follows: 

Volume  I Hardware  Design 

Volume  III  PAVSS  Operation  and  Software  Documentation 
Volume  IV  PAVSS  Program  Summary  and  Recommendations 


PREFACE 


The  problems  related  to  aircraft  trailing  vortices  are  currently  under 
intensive  study  for  the  Federal  Aviation  Administration  (FAA)  by  the 
U.  S.  Department  of  Transportation.  The  Transportation  Systems  Center 
(TSC)  of  the  U.  S.  Department  of  Transportation  (DOT)  initiated  and  is 
carrying  out  several  programs  in  this  area,  including  programs  to 
develop  acoustic  systems  for  detecting,  tracking,  and  measuring  the 
strength  of  aircraft  wake  vortices.  Avco  Corporation's  Systems 
Division  (Avco/SD)  designed,  built,  and  tested  a pulsed  acoustic 
vortex  sensing  system  (PAVSS)  under  Contract  DOT-T&C---620  in 
support  of  these  DOT /TSC  efforts.  The  software  currently  used  in 
vortex  tracking  produces  useful,  but  difficult-to-interpret  data  tracks. 
Avco/SD  engaged  two  subcontractors  to  investigate  possible  improve- 
ments in  the  area  of  vortex  tracking  using  pulsed  acoustic  data.  This 
report  presents  Avco's  recommendations  regarding  incorporation  of 
various  improvements  proposed  by  the  subcontractors,  and  includes 
the  final  reports  prepared  by  each  of  the  subcontractors. 

This  volume  of  the  final  report  on  the  PAVSS  program  describes  the 
approaches  made  toward  improving  the  software  used  for  vortex 
tracking-  Other  aspects  of  the  system  are  covered  in  additional  volumes 
(Volumes  I,  III,  and  IV). 

The  work  performed  under  this  contract  was  significantly 
enhanced  by  the  close  cooperation  and  contributions  of  Ralph  Kodis, 

David  Burnham,  and  Thomas  Sullivan,  all  of  DOT/TSC. 
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INTRODUCTION 


This  introduction  is  presented  in  two  parts.  Paragraph  1. 1 serves  as 
an  introduction  to  the  complete  final  report;  Paragraph  1.2  serves  as 
an  introduction  to  this  volume  (Volume  II)  of  that  report. 


1.  1 INTRODUCTION  TO  THE  FINAL  REPORT 


Trailing  vortices  from  heavy  jet  aircraft  represent  a currently  undefined 
hazard,  particularly  during  landing  and  takeoff  operations.  Considerations 
of  safety  and  the  need  to  optimize  airport  operation  make  it  essential  to 
acquire  positive  information  about  the  presence  and  locations  of  vortices 
generated  by  heavy  aircraft. 

The  feasibility  of  using  multi-static  pulsed  acoustic  radar  to  detect  and 
track  wake  vortices  has  been  demonstrated  by  the  Department  of 
Transportation's  Transportation  Systems  Center  (DOT/TSC)  in  tests  at 
Logan  International  Airport,  Boston,  Mass.;  at  John  F.  Kennedy  Inter- 
national Airport,  New  York,  N.  Y.;  and  at  the  National  Aviation  Facilities 
Experimental  Center  (NAFEC),  Atlantic  City,  N.  J.  The  hardware  used 
during  these  tests  consisted  of  laboratory  models.  The  equipment  was 
not  engineered  for  long-term  installation  in  the  field,  and  was  incapable 
of  automatic  real-time  data  processing  and  display. 


This  report  describes  a pulsed  acoustic  vortex  sensing  system  (PAVSS) 
development  program  carried  out  by  Avco  Systems  Division  (Avco/SD) 
for  DOT/TSC  under  Contract  DOT-TSC-620.  The  goal  of  this  program 
was  to  develop,  build,  and  test  an  engineered  wake  vortex  sensing  system 
consisting  of  acoustic  sensors  and  associated  electronics;  to  acquire  and 
process  the  sensed  data;  and  to  display  this  data  visually  in  real  time. 


The  complete  final  report  on  this  program  consists  of  this  volume 
(Volume  II,  STUDIES  OF  IMPROVED  PAVSS  PROCESSING  TECHNIQUES) 
and  three  additional  volumes,  as  follows; 

Volume  I HARDWARE  DESIGN 


Volume  III  PAVSS  OPERATION  AND  SOFTWARE 

DOCUMENTATION 

Volume  IV  PAVSS  PROGRAM  SUMMARY  AND 

RECOMMENDATIONS 
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1.2  INTRODUCTION  TO  VOLUME  U,  STUDIES  OF  IMPROVED 
PAVSS  PROCESSING  TECHNIQUES 

This  volume  of  the  Final  Report  on  the  Pulsed  Acoustic  Vortex  Sensing 
System  describes  two  studies  carried  out  in  support  of  that  program. 
These  were  independent  studies,  carried  out  by  Scope  Electronics,  Inc. 
and  the  Arcon  Corporation  under  Avco/SD  subcontracts,  to  investigate 
methods  of  improving  automatic  vortex  delay  tracking  techniques.  The 
final  reports  by  these  subcontractors  are  included  in  this  volume  as 
Appendixes  A and  B,  respectively. 

a.  Background 

Avco/SD  designed  and  developed  a PAVSS  capable  of  detecting 
and  locating  trailing  (wake)  vortices  produced  by  aircraft 
during  landings.  In  achieving  this  capability,  the  system  must 
track  a series  of  acoustic  pulses  that  have  been  refracted  by  a 
vortex  (the  vortex  return).  Such  returns  are  slightly  delayed 
relative  to  direct  (line-of-sight,  LOS)  acoustic  pulses.  Avco/ 

SD  developed  a simple  tracking  algorithm  designed  to  track 
both  the  vortex  returns  and  the  LOS  returns  in  the  presence 
of  noise.  This  approach,  while  meeting  the  goal  of  simplicity, 
presented  certain  problems  related  to  initial  acquisition  of  a 
vortex  return  track,  to  tracking  a vortex  return  close  (in  time) 
to  an  LOS  signal,  to  bridging  gaps  in  a vortex  return  track,  and 
to  stopping  the  track  upon  reaching  the  end  of  vortex  return  data. 

Avco/SD,  therefore,  placed  subcontracts  with  two  companies 
(Scope  Electronics,  Inc.  , of  Reston,  Virginia,  and  Arcon 
Corporation,  of  Wakefield,  Mass.,  Avco/SD  Purchase 
Orders  259700  and  259699,  respectively)  to  conduct  inde- 
pendent investigations  of  tracker  development.  This  was 
done  to  take  full  advantage  of  the  current  state-of-the-ait 
expertise  of  these  firms  in  the  area  of  tracker  technology. 

Appendixes  A and  B are  the  final  reports  by  Scope  (Report 
SEI-7063,  dated  10  May  1974)  and  Arcon  (Report  R74-2W, 
dated  18  March  1974). 

b.  Organization  of  Volume  II 

The  remaining  sections  of  this  volume  of  the  final  report 
cover  the  areas  described  below: 

Section  2 Discusses  the  major  features  of  the 

investigation  and  approach  followed 
by  each  subcontractor. 
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Section  3 Presents  Avco/SD's  recommendations 

regarding  implementation  of  certain 
tracker  features. 

Appendix  A Final  Report  by  Scope  Electronics,  Inc 

Appendix  B Final  Report  by  Arcon  Corporation. 
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2.  DISCUSSION 


This  section  summarizes  the  features  of  the  approach  that  each  sub- 
contractor employed.  Details  regarding  their  efforts  are  contained 
in  the  appropriate  appendixes  (Appendix  A for  Scope,  Appendix  B for 
Arcon  which  are  presented  essentially  as  received). 

2.  1 SCOPE  ELECTRONICS,  INC. 

Scope  Electronics,  Inc.  started  with  the  basic  ‘'minimum  mean 
square  error"  tracker  developed  by  Avco/SD  for  the  PAVSS-  program 
and  improved  it  in  certain  specific  areas.  The  modifications  consist 
essentially  of  a special  starting  routine,  tighter  tracking  with  a non- 
symmetric  tracking  window,  line-of-sight  masking,  and  a temporary 
stop  routine.  These  improvements  were  implemented  and  tested 
(using  taped  data  covering  15  different  landings).  The  results  indicated 
that  the  modified  program  showed  significant  improvement  in  each  of 
the  problem  areas. 

The  principal  features  of  the  Scope  program  are; 

1.  Tracking  of  the  LOS  data  prior  to  aircraft  arrival  with 
appropriate  guard  bands. 

2.  Use  of  a delay  tracker  with: 

a.  A start  routine  that  employs  a bin  approach. 

b.  Side  track  avoidance  capability. 

c.  Provisions  for  re-starting  when  track  gaps  occur. 

d.  An  effective  stop  routine. 

Except  for  the  LOS  guard  band  provisions,  all  of  the  features  proposed 
by  Scope  are  new  and  are  directly  applicable  to  the  PAVSS  minicomputer 
software. 

2.  2 ARCON  CORPORATION 

The  Arcon  Corporation  developed  an  entirely  new  tracking  technique, 
one  based  upon  use  of  a fixed-point  smoothing  algorithm.  The  algorithm 
represents  an  implementation  of  a Kalman-Bucy  filter  based  upon  an 
idealized  model  of  track  dynamics.  Due  to  the  fixed-point  smoothing 
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algorithm's  efficiency,  multiple  tracking  is  possible,  and  its  use  over- 
comes the  various  tracking  difficulties  that  were  encountered  using 
Avco's  single-track  technique.  A simplified  version  of  the  Arcon 
method  has  been  programmed  and  checked  out.  The  program  has  only 
been  given  preliminary  tests  (with  promising  results).  No  formal  test 
results  were  obtained. 

The  Arcon  Corporation  program  features: 

1.  Multiple  tracking,  using: 

a.  Parent  and  trial  tracks. 

b.  Track  selection  based  on  vortex  calculations. 

2.  A smoothing  approach  based  on  the  Kalman-Bucy  filter. 

Since  Arcon  provided  no  test  runs  using  these  features,  their  actual  value 
in  alleviating  the  present  PAVSS  software  difficulties  could  not  be 
determined. 
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Avco/SD  recommends  that  consideration  be  given  to  implementing 
several  of  the  tracker  improvements  developed  in  the  course  of  the 
Scope  and  Arcon  investigations. 

Initial  test  results  obtained  with  the  present  Avco/SD-built  PAVSS 
and  its  associated  software  indicate  that  the  most  pressing  needs  are 
for;  (1)  improved  start,  or  delay  track  acquisition,  a~d  (2)  better 
smoothing  of  the  track  after  it  is  acquired.  Track  smoothing  is  also 
important  from  the  viewpoint  of  production  of  more  meaningful  and 
more  easily  interpreted  plot  presentations.  The  features  incorporated 
in  the  program  proposed  by  Scope  appear  to  satisfy  all  of  these  needs. 

The  bin  approach  to  the  starting  routine,  which  is  initiated  when  air- 
craft noise  has  subsided,  should  provide  excellent  initial  acquisition 
of  the  delay  track.  The  side-track  avoidance,  re-start,  and  stop 
routines  should  insure  that  the  best  delay  track  is  maintained  ' nile 
the  vortex  remains  in  the  receiver /transmitter  field  of  view.  This 
should  prevent  acquisition  of  ground  clutter  during  acquisition  of 
weak  or  intermittent  delay  track  data. 

Avco/SD  strongly  recommends  that  the  PAVSS  software  be  modified 
immediately  to  implement  the  features  discussed  above,  thereby 
substantially  improving  the  system's  tracking  capability. 

The  Kalman- Bucy  filter  approach  investigated  by  Arcon  had  pre- 
viously been  considered  by  both  Avco/SD  and  TSC.  No  test  data 
has  been  supplied  to  evaluate  the  Arcon  technique.  It  may,  however, 
prove  desirable  to  employ  it  in  several  test  cases  to  permit  its 
evaluation. 
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Figure  1.  System  Configuration 


Section  1 


INTRODUCTION 


This  is  the  final  report  on  the  automatic  detection  and  tracking 
of  the  pulsed  acoustic  vortex  data  after  the  receiver  processing 
and  analog  to  digital  (A/D)  conversion.  The  details  of  the 
analog  signal  processing  and  A/D  conversion  prior  to  the  auto- 
matic tracking  of  the  vortex  from  the  digital  data  are  dealt 
with  in  the  AVCO  Pulsed  Acoustic  Vortex  Sensing  System  (PAVSS)  De- 
sign Evaluation  Report  (1).  Bistatic  Acoustic  Radar  Principles 
and  the  experimental  results  on  which  the  PAVSS  is  based  are  well 
documented  in  the  DOT-TSC  Technical  Reports  of  6/1971,  1/1972  and 
12/1972,  Ref.  (2),  (3),  (4)  respectively. 


For  the  sake  of  completeness,  a brief  description  of  the  PAVSS 
and  the  front-end  receiver  processing  will  be  given  where  applicable 
to  vortex  tracking.  A typical  configuration  of  the  two  vortices 
generated  by  an  aircraft  and  the  six  transceivers  used  to  transmit 
the  acoustic  pulse  and  receive  the  direct  line  of  sight  (LOS) 
return  and  the  scattered  delayed  return  is  shown  in  Figure  1. 

Vortex  Vj  will  scatter  the  pulse  from  TR1 , TR2,  and  TR3  forwards 
and  vortex  V2  will  scatter  the  pulse  from  TR4 , TRS,  and  TR6  forwards. 
Figure  2 shows  a typical  return  for  an  ideal  (noise-free)  receiver 
TR4.  The  analog  processing  of  the  receiver  output  and  A/D  conversion 


1 


are  summarised  in  the  following  steps.  For  details,  the  rcrder 

is  referred  to  Ref.  (1).  pp.  4-35  to  4-44. 

1)  High-pass  prefilter  the  signal  to  remove  the  low-frequency 
noise  without  attenuating  the  signal  component. 

2)  Pass  through  a synchronous  detector  which  acts  as  both  a 
narrow-band  filter  and  a detector. 

3)  Find  true  derivative  and  reject  negative  derivative. 

4)  Find  true  integration  (of  positive  derivative  only)  and 
reset  to  zero  each  time  negative  derivative  is  encountered. 
(The  overall  effect  of  3 and  4 is  to  restore  the  signal  but 
without  the  portions  with  negative  slopes  and  start  all 
portions  with  positive  slopes  from  zero  baseline.) 

5)  Compare  with  preselected  thresholds  and  mark  1 (dot)  for 
crossings  with  positive  slopes.  In  the  ideal  receiver, 
the  first  crossing  would  be  due  to  the  LOS  and  the  second 
due  to  the  delayed  return  from  the  closest  transmitter. 

(The  preselected  threshold  is  chosen  on  the  basis  of  noise 
statistics,  by  computing  a weighted  sum  of  the  noise  mean 
and  noise  variation  as  in  (5)  pp.  3-79.) 


Steps  1 through  5 are  shown  diagrammatically  in  Figure  3. 

This  processed  and  digitized  receiver  data  forms  the  basis  for 
automatic  detection  and  tracking  of  the  line  of  sight  (LOS)  and 
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vortex  returns.  Knowing  the  geometry  of  the  transceivers  and 
the  speed  of  sound,  the  LOS  and  delay  tracker  results  are 
then  used  to  compute  the  vortex  position.  Redundancy  is  built 
into  the  system  because  under  ideal  conditions  only  2 receivers 
(transmitters)  and  1 transmitter  (receiver)  are  required  to  locate 
the  vortex  position  which  is  given  by  the  intersection  of  two 
ellipses  with  one  common  focus  at  the  single  transmitter  (receiver) 
and  the  other  foci  at  the  two  receivers  (transmitters)  as  shown 
in  Appendix  C and  D in  Ref.  (5) . A system  diagram  is  shown  in 
Figure  4.  This  effort  was  restricted  only  to  the  part  block 
diagrammed  as  "DELAY  TRACKER."  Magnetic  tape  consisting  of  pro- 
cessed and  digitized  data  from  15  aircraft  landings  were  provided 
for  testing  the  performance  of  the  "DELAY  TRACKER". 

A few  points  regarding  the  analog  to  digital  (A/D)  conversion 
should  be  noted.  The  amplitude  and  pulse  width  information  of 
the  analog  return  is  lost,  the  threshold  level  is  a constant  and 
time  sampling  precludes  getting  more  than  one  dot  (1)  in  a sampling 
interval  (approximately  .5  msec.)  in  a given  frame  for  all  receivers. 
A typical  receiver  output  after  A/D  conversion  is  displayed  in 
Figure  5 in  the  form  of  a delay  (approximately  .5  msec,  increments) 
vs.  time  (approximately  .5  sec.  increments)  plot. 
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Section  2 
THE  TRACKER 


2.1  The  Approach 

The  receiver  output  and  the  analog  processing  done  on  it  before 
digitization  were  briefly  discussed  in  the  previous  section. 

A typical  receiver  output  after  digitization  is  shown  in  Figure  5. 
as  a plot  of  delay  (approximately  .5  msec.)  vs.  time  (approximately 
.5  sec . ) . 

Our  approach  to  tracking  was  to  divide  the  output  into  four  dis- 
tinct non-overlapping  regions  functionally  labeled  as  (1)  LOS 
Track  Only,  (2)  Noise,  (3)  Start  and,  (4)  Delay  Track  Only. 

These  regions  shown  in  Figure  6 are  categorized  by  the  type  of 
tracking  done  in  them,  i.e.,  LOS  tracking  only  before  aircraft 
arrival  in  (1),  no  tracking  in  aircraft  noise,  (2)  Start  routine 
only  for  the  vortex  delay  track  implemented  in  (3) , and  Vortex 
delay  tracking  only  in  (4).  These  forms  of  tracking  will  be 
discussed  below. 

2.2  LOS  Tracking  Only 

Before  the  aircraft  arrival,  the  tracks  are  relatively  noise-free 
and  the  regions  where  the  LOS  tracks  are  expected  can  be  closely 
specified.  By  restricting  the  LOS  tracking  to  these  tight  regions, 
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the  problem  of  tracking  clutter  can  be  easily  overcome.  For  each 
LOS  track  there  are  two  quantities  of  interest  which  must  be 
extracted  for  later  vortex  tracking  and  calculations,  i.e.,  the 
LOS  running  average  and  LOS  running  variation.  Throughout  this 
discussion  the  term  variation  will  be  used  to  describe  a quantity 
similar  to  variance  but  differing  in  that  it  uses  the  absolute 
value  of  the  error  rather  than  the  square  of  the  error  as  shown 
in  (2.2).  This  contributes  to  great  savings  in  computation  with- 
out degrading  the  performance. 


Since  the  running  average  with  a time  constant  of  M would  require 
the  storage  of  M values,  in  general,  a slightly  different  approach 
yielding  a similar  answer  but  requiring  one  storage  element  was 
found  to  be  more  appropriate  as  shown  in  (2.1)  and  (2.2).  Subscript 
i represents  the  ith  frame  and  Y^  the  LOS  value  in  ith  frame. 


LOS  Running  Average  YS^: 


YSf 


M . YS^  + Yi 

ETT 


LOS  Running  Variation  YSS^: 

M . YSSi_1  ♦ |YSi  - Y. | 


YSSi 


M ♦ 1 


(2.1) 


(2.2) 


A-15 


10 


The  LOS  tracking  is  terminated  as  soon  as  the  noise  count  in 
each  frame  exceeds  a preset  threshold  based  on  the  difference 
between  aircraft  noise  count  and  the  count  in  the  relatively 
noise-free  region  prior  to  the  aircraft  arrival.  The  values  of 
YS^  and  YSS^  are  accepted  as  the  final  values  for  the  LOS  and 
LOS  variation  respectively  for  use  in  the  vortex  delay  tracking 
and  other  vortex  related  calculations.  As  for  the  aircraft  noise 
region  ideally  no  tracking  at  all  is  done  during  this  interval. 

Experimentally  it  was  found  that  a region  of  20  increments 
(approximately  10  msec.)  on  either  side  of  the  expected  LOS  was 
more  than  sufficient  for  tracking  LOS.  Also  the  time  constant 
M ■ 40  was  implemented.  It  was  further  found  that  the  results  are 
not  degraded  even  if  the  tracking  is  continued  into  the  aircraft 
noise  region.  This  is  true  because  of  the  fact  that  the  time 
constant  M » 40  is  much  larger  than  the  aircraft  noise  region. 

As  far  as  the  experimental  results  are  concerned  this  was  the 
procedure  that  was  followed  with  the  above  parameters.  They  are 
shown  in  Figures  17  through  31  with  suffix  b as 

NF  « - LOS  «=  - LOSV  = - 

where 

NF  = Total  Number  of  Frames  used  for  averaging 

LOS=  Final  LOS  average  in  approximately  .5  msec  increments 

LOSV  ■ Final  LOS  Variation  in  approximately  .5  msec,  increments 
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2.3  Running  MMSE  Delay  Tracker 


Once  the  aircraft  noise  has  subsided  significantly  (ascertained 
by  the  noise  count  in  a frame)  it  is  possible  to  do  the  delay 
tracking.  Since  the  LOS  mean  and  variation  were  finalized  prior 
to  the  aircraft  arrival  they  will  be  assumed  to  remain  unchanged 
during  delay  tracking.  This  approach  not  only  reduces  the  num- 
ber of  computations  required  during  delay  tracking,  (essential 
for  real-time  tracking  with  a minicomputer),  but  also  keeps  the 
LOS  mean  and  variation  free  of  any  vortex  interference.  This  is 
especially  true  for  low-lying  vortices  close  to  the  LOS. 

Th»  quantities  required  for  the  running  minimum  mean  square  error 
(MMSE)  delay  tracker  are  shown  in  Table  1 for  clarity.  The 
computations  required  for  the  MMSE  delay  tracker  are  tabulated 
in  Appendix  A.  Again,  as  in  the  LOS  calculations,  it  is  convenient 
to  deviate  slightly  from  the  exact  running  average  in  order  to 
save  on  the  storage  requirements.  The  exact  running  average  would 
require  a storage  of  as  many  quantities  as  the  time  constant 
(e.g.  N in  the  case  of  c^) , however,  the  modified  approach  requires 
only  one  storage  element  per  running  average  as  shown  in  Appendix 
A. 
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RCOUSTOGRflM  CH.  5 


The  basic  delay  tracker  as  implemented  by  AVCO  was  found  to 
suffer  from  the  5 problems  categorized  below: 


1)  Start 

2)  Side-track 

3)  LOS  Masking 

4)  Clutter 

5)  Stop 

2. 3. 1.1  Start 

Start  refers  to  the  difficulty  in  acquiring  a delay  track  at  the 
beginning.  The  twin  problem  of  acquiring  the  right  track  and 
acquiring  it  as  fast  as  possible  was  found  to  be  the  most  important 
flaw  in  the  basic  delay  tracker.  An  example  of  the  start  problem, 
is  shown  in  Figure  7. 

The  solution  to  the  start  problem  was  to  divide  the  vortex  delay 
into  two  regions,  the  first  region  would  be  exclusively  used  for 
implementing  the  Start  Routine  and  the  following  region  for  delay 
tracking.  The  Start  Routine  was  to  accomplish  the  delay  track 
acquisition  correctly  and  as  fast  as  possible  with  no  more  com- 
putations required  than  those  used  for  the  delay  tracking.  To 
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this  end  the  idea  selected  was  to  divide  the  receiver  region 
above  the  LOS  into  horizontal,  sloping  and  offset  sloping  bins 
as  shown  in  Figure  8. 

A score  of  the  number  of  hits  occurring  in  each  bin  is  kept 
until  one  of  them  crosses  a preselected  threshold.  In  case  of 
a draw  between  the  horizontal  bin  and  the  sloping  bins  a prefer- 
ence is  given  to  the  sloping  bins.  When  two  or  more  bins  with 
the  same  slope  tie,  the  uppermost  bin  is  selected. 

The  delay  point  within  the  winning  bin  and  its  slope  are  then 
used  to  initiate  the  delay  tracking.  All  the  initial  values  of 
the  tracking  parameters  in  Appendix  A are  computed  as  if  the 
number  of  hits  in  the  winning  bin  occurred  consecutively  and 
with  the  slope  of  the  bin  as  shown  in  Appendix  B.  This  assures 
continuity  and  proper  weighting  of  the  past  delay  data. 

As  far  as  £he  implementation  of  the  Start  Routine  is  concerned, 
a bin  width  of  10  increments  or  about  5 msec  was  found  to  be 
optimum.  It  was  found  convenient  to  ignore  the  contents  of  the 
first  2 horizontal  bins  and  4 sloping  and  offset  sloping  bins 
immediately  above  the  LOS  because  of  the  concentration  of  noise 
close  to  the  LOS  and  the  fact  that  the  vortices  generally  require 
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reasonable  amount  of  time  before  they  descend  to  such  low  heights. 
Because  of  the  ease  with  which  vortices  with  large  delays  were 
acquired  and  their  rapidly  descending  nature  and  the  available 
computer  storage  it  was  found  sufficient  to  include  only  8 
horizontal  and  16  sloping  and  16  offset  sloping  bins.  Also  the 
Start  Routine  was  terminated  if  the  threshold  was  not  crossed 
after  searching  through  60  frames  (approximately  30  sec). 

Results  of  the  above  implementation  are  presented  in  Figures  17 
through  21. 

All  of  these  parameters  may  be  changed  if  the  situation 
warrants  it.  If  the  storage  requirements  of  the  minicomputer 
necessitate  it,  the  elimination  of  Offset  Sloping  Bins  would 
not  affect  track  acquisition  significantly,  however,  it  would 
affect  the  total  time  required  for  track  acquisition. 


2.3. 1.2  Side-track 

Side-tracking  refers  to  the  problem  of  getting  off  the  main 
delay  track  and  tracking  a ghost  image  as  shown  in  Figure  9 or 
tracking  into  a false  track  composed  of  noise  as  shown  in  Figure 
10.  This  difficulty  arises  due  to  the  fact  that  when  a hit 
belonging  to  the  correct  delay  track  does  not  occur,  the  tracker 
picks  any  other  closest  hit  which  occurs  within  the  frame  and 
starts  tracking  along  the  new  path. 


SRMPLE 

too  500  600  700  600 


SAMPLE 


The  solution  to  this  problem  was  to  reject  any  hits  which  occur 
outside  a preselected  window  around  the  predicted  estimate  of 
the  delay.  Since  the  Start  Routine  enables  us  to  pick  the  cor- 
rect track  with  greater  confidence,  we  can  afford  to  put  a 
tighter  bound  around  the  delay  tracker.  To  allow  the  tracking 
of  upward  rising  vortices  the  window  is  made  wider  on  the  upprr 
side  as  shown  in  Figure  11.  A variable  window  based  on  the 
delay  tracker  variation  was  also  tried  and  found  to  be  inappro- 
priate for  the  type  of  vortex  data  provided.  The  variable  win- 
dow which  has  a tendency  to  get  very  narrow  as  the  delay  tracker 
goes  through  smooth  data  quite  often  loses  the  track  during 
transition  from  smooth  to  rough  data  and  therefore  this  idea 
was  abandoned. 

As  far  as  the  implementation  is  concerned,  a window  of  7 incre- 
ments (approximately  3.5  msec)  on  either  side  with  an  extra 
margin  of  7 increments  on  the  upper  side  was  found  to  give  satis- 
factory results.  These  values  may  be  changed  slightly  without 
affecting  the  results  significantly;  however,  it  was  found 
necessary  to  retain  the  asymmetry  in  the  window  to  track  the 
upward  rising  vortices  without  substantially  affecting  the 
tracking  of  normal  descending  vortices.  Some  of  the  results 
can  be  seen  in  Figures  22  through  24. 
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2. 3. 1.3  LOS  Masking 

LOS  masking  refers  to  the  problem  of  tracking  the  vortex  delay 
as  close  to  the  line  of  sight  as  possible.  This  is  especially 
true  for  the  low-lying  vortices  close  to  ground  which  persist 
for  a long  period  of  time.  The  problem  of  delay  tracks  which 
interact  with  the  LOS  and  disappear  into  the  LOS  only  to  reappear 
later  as  distinct  delay  tracks  will  be  dealt  with  in  the  Stop 
Routine. 

The  solution, which  gave  favorable  resultsfwas  to  select  a region 
above  the  LOS  based  on  the  LOS  statistics  which  would  be  off 
limits  to  the  delay  tracker  as  shown  in  Figure  12.  The  region 
was  set  equal  to  two  times  the  LOS  variation  or  5 increments 
(approximately  2.5  msec.)  whichever  was  maximum.  The  minimum  value 
of  5 increments  appears  to  be  a reasonable  number  to  allow  for 
the  interaction  of  the  vortex  with  the  LOS. 

The  above  criterion  together  with  the  tighter  asymmetrical  window 
with  a smaller  value  on  the  lower  side,  keeps  the  delay  tracker 
from  getting  pulled  by  the  LOS  and  yet  there  is  enough  flexibility 
to  track  vortices  which  may  rise  again. 

In  the  implementation,  the  off  limits  region  above  the  LOS  was 
chosen  as  the  maximum  of  two  times  the  LOS  variation  or  5 incre- 
ments and  the  lower  window  was  set  to  7 increments  (approximately 
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3.5  msec).  The  results  of  the  LOS  masking  (some  with  the  STOP 
routine)  are  shown  in  Figures  25  through  28. 

2. 3. 1.4  Clutter 

Clutter  problems  can  be  divided  into  two  categories,  clutter 
within  the  expected  region  for  delay  tracking  and  clutter  out- 
side this  region.  In  the  15  landings  that  were  examined  no 
obvious  cases  of  the  former  were  noticed.  A suggested  solution 
for  the  clutter  within  the  delay  tracking  region  would  be  to 
track  these  before  the  aircraft  arrival  in  the  segment  marked 
LOS  Tracking  Only  in  Figure  6.  Once  the  clutter  tracks  are 
known  they  can  be  used  in  the  Start  Routine  by  eliminating  those 
horizontal  bins  which  contain  clutter  from  consideration.  Once 
the  correct  start  is  made,  the  tighter  window  around  the  delay 
tracker  will  keep  it  from  being  sidetracked. by  clutter. 

As  far  as  the  clutter  outside  the  expected  delay  tracking  region 
is  concerned,  these  can  conveniently  be  accomplished  by  masking 
out  those  regions  completely.  In  the  implementation  it  was 
found  that  the  horizontal  bins  above  120  increments  (approxi- 
mately 60  msec)  from  the  LOS  were  not  required,  hence  clutter 
in  this  region  would  not  pose  any  problems.  Figures  31  and  32 
show  the  results. 
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2 . 3 . 1 . 5 Stop 

Finding  suitable  criteria  for  stopping  the  tracker  was  another 
problem  encountered  in  implementing  the  basic  delay  tracker. 

The  three  criteria  shown  in  Figure  13  as  1)  Track  Density, 

2)  Noise  Density  and  3)  Delay  Error  were  examined  for  their 
applicability. 

1)  Track  Density  refers  to  the  running  average  of  density  of 
hits  within  the  delay  tracker  window.  It  is  shown  as 

in  Appendix  A. 

2)  Noise  Density  refers  to  the  running  average  of  the  density 
of  frames  with  more  than  one  hit  above  the  LOS  masking 
region  per  delay  track. 

3)  Delay  Error  Tefers  to  running  average  of  the  absolute  error 
between  the  calculated  estimate  and  the  actual  delay  which 
is  shown  as  D in  Appendix  A. 

Except  for  the  noise  only  situations,  the  track  density  param- 
eter or  a slight  variation  thereof  was  found  to  be  quite  suf- 
ficient for  stopping  the  delay  tracker.  Since  the  track  density 
parameter  is  restricted  to  a tight  window  around  the  delay 
tracker,  the  delay  error  criterion  was  found  to  be  redundant 
for  this  case. 

The  problem  of  distinguishing  noise  only  situations  from  noise 
plus  delay  track  situations  as  shown  in  Figure  14  proved  to  be 
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Noise  Plus  Track 


NOISE  & TRACK 
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more  difficult.  Since  the  noise  density  criterion  as  stated 
earlier  (#2  in  Figure  13)  would  give  similar  results  in  both 
these  cases  and  since  the  tracker  always  chooses  the  path  with 
the  minimum  mean  squared  error  through  either  of  these  cases, 
it  would  appear  that  none  of.  the  criteria  in  Figure  13  would 
be  able  to  distinguish  between  these  two  cases. 

It  seems  that  the  noise  density  criterion  needs  to  be  modified. 

One  promising  approach  would  be  to  consider  the  noise  density  in 
the  tight  window  around  the  delay  tracker  instead  of  the  larger 
delay  frame.  This  would  keep  the  noise  density  for  the  noise  only 
case  relatively  unchanged  but  drastically  reduce  the  noise  density 
in  the  noise  plus  delay  track  case. 

As  far  as  the  situation  where  the  delay  track  disappears  in  the 
LOS  and  reappears  later  to  hover  above  the  LOS  is  concerned,  it 
was  found  convenient  to  make  the  first  stop  a temporary  one  and 
initiate  a Restart  Routine.  The  Restart  Routine  consists  of  looking 
through  the  tighter  window  for  a preselected  number  frames  until  the 
number  of  hits  crosses  a certain  threshold  which  in  general  is  the 
same  as  that  in  the  Start  Routine.  If  the  number  of  hits  does 
not  cross  the  threshold  then  the  delay  tracker  is  terminated 
permanently.  If  the  number  of  hits  crosses  the  threshold,  then 
the  delay  tracking  is  resumed  until  the  second  stop  which  termi- 
nates the  delay  tracker  permanently. 
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A slight  modification  of  this  would  be  to  reinitiate  the  Start 
Routine  after  the  first  temporary  stop.  Since  the  temporary 
stop  generally  occurs  close  to  the  LOS  and  noise  is  also  more 
pronounced  close  to  the  LOS,  the  Start  Routine  may  initiate 
tracking  of  noise  only.  This  problem  may  be  alleviated  some- 
what by  1)  choosing  only  those  Start  Routine  bins  which  are 
above  the  lower  bound  of  the  temporarily  stopped  delay  tracker, 
2)  assigning  a different  slope  to  the  Sloping  Bins,  and  3)  de- 
creasing the  bin  width  since  fewer  bins  are  required. 

In  Figures  17  through  32  a slight  modification  of  track  density 
criterion  was  implemented.  The  delay  tracker  was  temporarily 
stopped  when  there  were  no  hits  in  8 consecutive  frames  within 
the  tracker  window.  The  tracker,  however,  continues  looking 
through  the  window  and  updating  all  the  parameters  and  if  the 
number  of  hits  within  the  window  exceeds  6,  restarts  tracking 
until  8 consecutive  no  hits  in  the  window.  Also  if  the  number 
of  hits  within  the  delay  tracking  window  does  not  exceed  the 
threshold  of  6 in  the  remaining  frames,  the  Restart  Routine 
stops  the  delay  tracker  permanently.  In  general,  a number  such 
as  60  consecutive  frames  (approximately  30  sec)  would  be  more 
reasonable  for  the  Restart  Routine  implementation.  The  results 
of  the  Stop  Routine  implementation  can  be  seen  vividly  in 
Figures  26  through  30. 


A-35 


t 


30 


_ | | — ™ — 

'"7\  ‘ 


2.3.2  Flow  Charts  and  Computational  Requirements 

Since  the  basic  tracker  is  the  same  as  that  implemented  by  AVCO 
(1),  i.e.,  minimum  mean  square  error  delay  tracker,  its  flow 
chart  computational  requirements  and  timing  will  not  be  dis- 
cussed here.  The  flow  chart,  time  and  storage  requirements  for 
all  the  additions  and  modifications  discussed  in  the  previous 
sections  will  be  the  subject  of  this  section. 

2. 3. 2.1  Start  Routine  Implementation 

The  flow  chart  for  the  Start  Routine  is  shown  in  Figure  15.  The 
first  conditional  branch  marked  "START=1"  decides  whether  the 
Start  Routine  or  the  delay  tracker  is  in  effect.  START  is  of 
course  initialized  to  0 beforehand.  The  box  computes  the  bins 
in  which  the  delay  point  lies,  i.e.,  J2  for  horizontal,  J4  for 
sloping  and  J6  for  offset  sloping,  and  then  adds  one  to  the  pre- 
vious number  of  hits  in  the  appropriate  bin,  i.e.,  HITSH  for 
horizontal,  HITSS  for  sloping  and  HITSO  for  offset  sloping.  The 
following  three  conditional  branches  update  the  maximum  number 
of  hits,  MAXH  if  any  of  the  latest  bins  exceed  the  previous  max- 
imum, otherwise  MAXH  remains  unchanged.  The  delay  and  the  slope 
of  the  bin  corresponding  to  the  maximum  is  also  updated  or  left 
alone  in  parallel  with  MAXH  as  shown  in  the  three  boxes  to  the 
right. 

At  the  end  of  the  frame  a conditional  branch  compares  the  maxi- 
mum number  of  hits  MAXH  thus  far  with  the  threshold  HTHR  to 
decide  whether  to  continue  Start  Routine  or  to  initiate  delay 
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tracker.  The  last  conditional  branch  compares  the  total  number 
of  frames  gone  through  so  faT  with  the  preselected  threshold  to 
decide  whether  to  continue  the  Start  Routine  or  Stop  permanently. 


The  time  and  storage  requirements  for  the  above  Start  Routine 
on  PDP-11  with  the  extended  Arithmetic  Unit  (KE11-A)  are  as 
follows . 


Time:  6 Conditional  branches 

4 Multiplies/Divides 
6 Adds 


approximately  45  usee,  per  track 


Storage:  8 Words  HITSH 

16  Words  HITSS 
16  Words  HITSO 
6 Words  Miscellaneous 


46  Words  total  per  track 


There  are  a few  observations  worth  noting  about  the  Start 

Routine. 

1)  It  is  inherently  fixed  point. 

2)  Although  the  storage  requirement  appears  to  be  high 
compared  to  the  storage  requirement  for  the  delay  tracker, 
the  time  requirement  is  much  lower. 

3)  Since  each  word  is  expected  to  be  less  than  4 bits,  these 
words  can  be  packed  4 to  1 in  the  16-bit  computer  word 
requiring  approximately  12  computer  words.  The  time  required 
for  scrambling  and  unscrambling  the  words  is  traded  for 
lower  storage. 

4)  Since  the  Start  Routine  and  delay  tracker  are  operating  in 
mutually  exclusive  intervals,  the  storage  can  overlay  that 
used  by  the  delay  tracker. 
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Figure  16.  Flow  Chart  for  Remaining  Modifications  to  Delay  Tracker 
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2.3. 2. 2 Delay  Tracker  Additions 


The  flow  chart  for  the  remaining  modifications  under  the 
categories  of  side-tracking,  LOS  Masking,  Clutter  and  Stop 
are  shown  in  Figure  16.  The  first  three-conditional  branches  com- 
bine the  solution  to  side-tracking  and  LOS  Masking.  The  first 
branch  eliminates  data  points  within  LOS  variation.  The  next 
two  branches  accept  points  within  the  tracker  window  around  the 
predicted  estimate,  The  following  two  conditional  branches 

on  the  left  choose  the  point  with  the  lowest  error,  dQl<j.  The 
following  right  branch  performs  all  the  delay  tracker  updating  when 
dQid  has  been  changed  from  its  initial  large  value,  otherwise  it 
goes  through  the  Stop  Routine.  The  following  conditional  branch 
compares  the  number  of  consecutive  no  hits  to  a threshold  and 
Stops  the  first  time  if  the  threshold  is  exceeded.  The  conditional 
branch  "STOP  = 1"  decides  whether  it  is  a first  temporary 
stop  or  second  final  stop.  STOP  is  of  course  initialized  to  0 
at  the  beginning. 

The  maximum  time  and  storage  requirements  for  the  above  delay 
tracker  modifications  on  PDP-11  with  the  Extended  Arithmetic 
Unit  (KE11-A)  are  as  follows: 

Time:  5 Conditional  Branches 

1 Subtract 

1 Absolute  Value 

2 Stores 

Approximately  23  nsec,  per  track 
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Storage:  LOSV 


4 words  per  track 


It  is  important  to  note  that  some  of  the  storage  shown  in 


the  flow  chart  i.e.,  low,  up,  are  common  to  all  tracks  and  that 


the  path  marked  "NO  HIT  CALCULATION"  requires  an  insignificant 
number  of  computations  compared  with  "HIT  CALCULATION"  and  also 


occurs  less  frequently.  Therefore,  the  above  time  and  storage 


requirements  represent  the  most  frequent  and  time  consuming 
computational  additions  to  the  basic  delay  tracker.  In  Ref.  (1) , 
the  time  required  for  the  basic  delay  tracker  calculations  was 
•stimated  as  300  ysec.  per  track  in  fixed  point  arithmetic,  so 
the  above  additions  which  are  also  in  fixed  point  arithmetic  re- 
present an  increase  of  approximately  81. 


2.3.3  Fixed  Point  Implementation  of  Tracker 


Thi*  section  deals  with  the  problems  involved  in  using  the  finite 
16-bit  word  length  of  the  minicomputers  such  as  PDP-11  to  do  the 
delay  tracker  computations.  It  should  again  be  emphasised  that 


all  the  additions  and  modifications  done  on  the  basic  tracker  are 


inherently  fixed  point  and  the  16-bit  word  sise  is  amply  sufficient 


for  accomplishing  the  Start  Routine  or  the  additions  to  the 
basic  delay  tracker.  Therefore,  we  shall  be  concern-  ’ 1 y ■ ifh 

the  basic  tracker  computations  in  fixed  point  arithmetic. 

It  can  be  safely  stated  that  except  for  the  two  second  order 
non-central  moments  in  Appendix  A,  most  of  the  remaining  com- 
putations can  be  done  in  fixed  point  arithmetic  by  scaling  the 
quantities  up  and  occasionally  rescaling  if  necessary.  Since  the 
16-bit  woid  allows  up  to  65536  levels  of  quantization,  this  should 
be  quite  sufficient  for  all  the  quantities  in  Appendix  A except 
XSSi  and  XYS^  as  shown  below: 


M_ 


XSS, 


XSSi.1  + i2 


Ms  + 1 


and, 


XYS. 


MSj-l  * 1 • Yi 

Ms  ♦ 1 


Maximum  values  of  i and  Y^  (normalized  to  LOS)  were  found  to 
go  as  high  up  as  300  frame  increments  and  120  delay  increments 


2 

respectively,  hence  i exceeds  the  16-bit  word  size  and  i . Y^ 
is  very  close  to  the  limit.  The  solution  to  this  problem  is 
to  use  centralized  second  order  moments  as  follows: 


Var{X)i 


Ms  . Var{X}i.1  + (i  - XS i)2 
Ms  + 1 


and, 


Cov{XY>i 


Ms  . Cov{XY}i.1  + (i  - XSi)(Yi  - YS^ 
Ms  + 1 


Maximum  values  of  (i  - XS^  and  (Y^  - YSA)  should  not  go  beyond 
the  time  constant  (=10  frames)  and  the  larger  tracker  window 
width  (=14  increments)  respectively.  Besides  keeping  the  values 
low,  it  also  reduces  the  number  of  multiplies  by  2 since  the 
slope  A is  directly  given  by  Cov/Var. 


As  far  as  the  track  density  running  average,  is  concerned 
since  it  ideally  represents  an  integral  number  of  hits  per  (N  + 1) 
frames,  it  should  not  require  any  more  quantization  than  (N  + 1) 
levels.  In  order  to  avoid  the  degradation  due  to  the  cumulative 
effect  of  division  every  change  of  frame,  the  following  procedure 
should  work: 


oi.1  - N 2 * = CNum)i-_l 


N ♦ 1 


N * 1 


I 


I 


(Num) i_1  - aj.j  + Bi 
N + 1 


(Num) i 
N + 1 


For  the  present  frame  i,  by  using  mainly  the  numerator  from  the 
previous  frame  calculation  the  deterioration  from  the  fixed  point 
division  is  avoided  in  the  iterative  process.  This  also  reduces 
the  number  of  multiplies  by  one  at  the  expense  of  increasing  the 
storage  by  one. 

Assuming  Mg  is  quantized  to  the  same  level  as  i.e.,  (N  + 1) 
or  approximately  4 bits,  and  XS^_j,  YSi_1  do  not  exceed  8 bit  word 
size  for  integral  number  of  frame  and  delay  increments,  that  still 
leaves  4 more  bits  for  finer  quantization  of  XS^  and  YS^. 

Although  the  fixed  point  version  cf  the  basic  tracker  was  not 
implemented  as  such  the  results  of  the  modified  delay  tracker 
still  hold  since  the  changes  applied  to  the  basic  tracker  are  all 
in  fixed  point  arithmetic. 

2.4  Results 

This  section  discusses  the  effect  of  the  various  modifications 
to  the  basic  delay  tracker  on  actual  field  data  after  analog 
processing  and  A/D  conversion  by  AVCO.  The  following  data  were 
examined. 
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No.  of  Landings  IS 

No.  of  Receivers/Landing  6 

No.  of  Delay  Tracks/Receiver  3 

Total  no.  of  Tracks  examined  270 

For  the  sake  of  brevity, of  the  270  delay  tracks  only  those  which 
gave  problems  with  the  basic  tracker  and  the  corresponding 
improvement  from  the  modifications  are  shown  in  Figures  17 
through  32.  From  the  figures  it  can  be  seen  that: 

1)  The  Start  Routine  with  the  Horizontal,  Sloping  and  Offset 
Sloping  Bins  significantly  improves  the  wrong  start  and 
late  start  problem. 

2)  Tighter  delay  tracking  eliminates  sharp  discontinuities, 
wrong  tracks  and  tracking  the  LOS. 

3)  Non-symmetr ic  window  around  delay  tracker  enables  tracking 
of  upward  rising  vortices  and  prevents  tracking  LOS. 

4)  Masking  of  a region  2 times  LOS  variation. allows  tracking 
as  close  to  LOS  as  statistically  possible. 

5)  Temporary  stop  allows  for  tracking  of  delay  tracks  which 
disappear  into  LOS  only  to  reappear  later  as  low-lying  tracks 
which  tend  to  persist. 

From  the  visual  observation  of  the  modified  delay  tracker  results 
of  the  270  tracks  it  can  be  stated  that  the  miss  rate  (or  tracking 
wrong  tracks)  was  2 to  3%.  An  example  of  miss  is  shown  in  Figure 
26b.  Most  of  the  misses  were  due  to  situations  similar  to  those 
shown  in  Figure  26b,  where  there  are  two  equally  strong  delay 
tracks.  Some  misses  occurred  due  to  wrong  start  through  noise. 
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Noisy  situations  where  there  is  no  visible  delay  track  or  the 
delay  track  is  buried  in  too  much  noise  accounted  for  12%  of 
the  total  270  cases.  Although  the  modified  noise  density 
criterion  was  net  implemented,  it  is  believed  that  its  use 
would  significantly  improve  the  rejection  rate  of  noise  only 
situations  without  worsening  the  miss  rate  of  2 to  3%.  Since 
there  is  a built-in  redundancy  in  the  PAVSS,  both  the  miss  rate 
and  the  false  alarm  rate  would  be  drastically  reduced  when  the 
redundant  tracks  are  omitted  from  consideration.  Theoretically, 
since  only  2 tracks  out  of  the  possible  18  are  nonredundant , 
the  above  percentages  could  be  reduced  by  a factor  of  as  much 
as  9,  assuming  uniform  distribution. 
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Improved  Track  fTcra  Figure  19a  (-»700  samples  into  frame) 
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Improved  Track  frcm  Figure  21a  ( — 100  samples  into  frame) 
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Figure  25a 


Improved  Track  from  Figure  25a  (—700  samples  into  frame) 
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Improved  Track  from  Figure  28a  (-—100  samples  into  frame) 
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Spurious  Response,  Figure  31a  (~230  samples  into  frame) 


Section  3 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  objectives  of  this  effort  were  to  come  up  with  an  algorithm 
for  the  tracking  of  LOS  and  vortex  delay  which  did  not  have  any 
of  the  problems  associated  with  the  basic  minimum  mean  square 
error  (MMSE)  tracker,  so  that  it  can  eventually  be  implemented 
in  real-time  on  a minicomputer  such  as  PDP-11.  In  spite  of  the 
fact  that  the  processed  and  digitized  data  did  not  contain  the 
amplitude  and  bandwidth  (or  pulse  width)  information,  which  were 
considered  to  be  serious  handicaps,  all  the  goals  of  this  effort 
have  been  accomplished  and  suggestions  for  improvement  and  further 
work  are  also  made. 

Starting  with  the  basic  MMSE  tracker,  all  the  problems  associated 
with  it  were  individually  analyzed  and  various  solutions  tried. 
Successful  solutions  consisting  of  the  Start  Routine,  tighter 
tracking  with  non-symmetric  window,  LOS  masking  and  temporary 
stop  were  implemented  and  tried  on  a wide  cross-section  of  track 
data  from  15  landings.  Furthermore  the  solutions  consisted  of 
simple  procedures  done  in  fixed  point  arithmetic  which  can  easily 
be  added  to  the  basic  tracker  with  realistic  storage  and  time  re- 
quirements. 


JZ.it.-  J:  - i . 
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For  further  improvements  in  the  rejection  rate  for  the  noise 
only  situations  it  is  recommended  that  the  modified  noise 
density  criterion  mentioned  in  Section  2. 3. 1.5  be  implemented. 

It  would  also  be  interesting  to  further  pursue  the  fixed  point 
arithmetic  suggested  in  Section  2.3.3  for  the  basic  tracker  by 
implementing  the  total  modified  tracker  in  fixed  point  arithmetic 
and  trying  it  on  the  same  15  landings  or  any  other  later  improved 
data.  Real-time  implementation  of  the  final  algorithm  coded  in 
machine  language  on  a minicomputer  such  as  PDP-11  would  also 
be  of  great  value. 
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Appendix  A 

RUNNING  MMSE  DELAY  TRACKER  CALCULATIONS 


N ■ Time  constant  of  running  average  of  track  density 
M * Time  constant  of  running  MMSE  tracker 
i » Frame  number 

Value  of  the  delay  point  in  frame  i 
■1  if  a hit  occurs  in  frame  i 

Bi> 

A "0  otherwise 


TRACK  DENSITY, 


N . Oj.j  + Bj 
N + 1 


TC  MODIFICATION, 


M#  ■ • . M 


FRAME  NUMBER, 


M,  « XSi-l  ^ 1 

M,  ♦ ! 


DELAY, 


M,  • YSi-l  * Yi 
M + 1 

3 


(FRAME) • (FRAME) 


Ms  . jMlml  + r 
Ms  ♦ 1 


" 


(FRAME) • (DELAY) , 


SLOPE, 


INTERCEPT, 


DELAY  ESTIMATE , 


DELAY  ERROR, 


XYSi-l  + i • Yi 
Ms  + i 


XS.  . YS. 
i i_ 

XSt  . XSA 


B = YS£  - XS£  . A 


Yi  = B + A . i 


Ms-Di-1  + lYi  - Yil 


Ms  + i 


Di  " Di-1 


, if  Ba  - 1 


, if  Bi  = 0 


Appendix  B 

INITIAL  VALUES  OF  THE  TRACKER  PARAMETERS 


Assume  the  Start  Routine  gives, 

1)  hit  threshold  ■ MAXH 

2)  winning  delay  - MAXI 

3)  winning  bin  slope  * SLOPE 

Track  Density,  aJN  ■ 1 

Frame  Number,  XSjjj  ■ (MAXH  + l)/2 

Delay,  YStm  - (MAXI  - SLOPE  . (MAXH  - l)/2 

I 


(Frame)  • (Frame),  XSSJN  - (MAXH+1) (2 *MAXH+l)/6 
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VORTEX  TRACKER  REPORT 


General  Status 

The  following  brief  report  is  an  account  of  a new  design  for  the 
delay  tracker  in  the  AVCO  pulsed  vortex  tracking  system.  The  phil- 
osophy of  this  design  is  to  produce  an  integrated  and  efficient  algo- 
rithm which  remedies  some  of  the  observed  performance  deficiencies 
of  the  current  system.  Among  these  difficulties  are: 

a.  Excessive  running  time  due  to  use  of 
floating  point  arithmetic. 

b.  Erratic  start-up  of  tracks. 

c.  Failure  to  correlate  data  with  track  in 
ambiguous  situations. 

The  new  system  is  a more  consistent  and  powerful  practical  method 
for  delay  tracking. 

Additionally,  a simplified  version  of  this  general  design  has  been 
programmed  and  checked  out  for  use  with  the  AVCO  delay  tracking 
simulator.  The  simplification  has  been  to  omit  the  secondary  search 
procedure  and  the  generation  of  split  tracks  in  ambiguous  data  situ- 
ations. The  program  has  undergone  preliminary  tests  (only  as  a means 
of  checkout)  and  results  are  encouraging.  No  formal  test  results  arc 
available  at  this  time. 

Further  development  of  this  program  would  require  parameter 
tuning  to  optimally  match  data  characteristics.  The  full  algorithm 
should  also  be  implemented  to  assess  the  effectiveness  of  its  additional 
features.  Extensive  testing  should  be  performed  to  validate  the  algo- 
rithm performance  with  real  data. 

Design  Outline 

The  following  design  description  is  organized  around  Figures  1 
through  13.  The  remaining  figures  provide  some  background  information. 
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Figure  I.  Delay  Tracker  Features 

All  fixed  point  processing  is  used  to  minimize  computer  running 
time.  The  smoothing  algorithm  is  such  that  special  scaling  to  pre- 
serve accuracy  is  unnecessary  and  it  may  be  programmed  in  fixed 
point  directly  as  written.  (Delays  from  zero  to  maximum  are  accommo- 
dated with  sufficient  accuracy  by  the  fixed  point  fractions  0 to  1.  ) 

The  smoothing  algorithm  is  computationally  very  efficient  and  esti- 
mates track  position  and  velocity  recursively.  It  is  a Kalman-Bucy 
filter  based  on  an  idealized  model  of  track  dynamics  but  should  b;  em- 
pirically useful  in  the  present  context. 

Each  receiver /transmitter  combination  is  handled  separat'  -/  for 
tracking.  For  each  combination  several  tracks  may  be  held  simulta1  - 
eously.  This  permits  tracking  of  ghosts  or  baseline  in  addition  to  the 
main  track;  but,  more  importantly,  if  one  track  inadvertently  follows 
false  data,  a second  track  can  simultaneously  acquire  and  follow  the 
true  data.  The  choice  between  them  is  made  on  the  higher  level  when 
their  data  are  selected  for  vortex  calculations. 

In  order  to  select  the  best  from  among  several  tracks,  running 
measures  of  the  status  and  quality  of  a track  must  be  kept.  These 
measures  are: 

a.  Ty^e  status  - The  track  is  being  initiated 
upon  (still  tentative),  it  is  undergoing  normal 
tracking  (normal),  or  it  is  a split  (parent  or 
trial). 

b.  Firmness  - A measure  of  the  past  consistency 
of  data-track  correlations  which  controls 
smoothing  constant  and  bin  size  selection. 

c.  Hit/Miss  Count  - A count  of  consecutive  hits 
for  initiation  decisions  or  of  consecutive  misses 
for  drop  decisions. 

Multiple  tracking  and  a variety  of  data  situations  are  handled  by 
an  extended  correlation  logic.  This  logic  implements  a continuing 
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auto-initiation  of  new  tracks  on  any  data  frame.  It  also  selects  the 
best  data  for  each  track  and  performs  all  track  status  decisions.  In 
particular,  when  data  of  a normal  track  is  missing  in  the  primary 
bin,  it  constructs  a larger  secondary  bin.  A unique  datum  in  the 
secondary  bin  leads  to  a split  trial  track,  while  the  original  track  is 
also  extrapolated  (the  parent).  This  situation  is  resolved  on  subse- 
quent frames. 


Figure  2.  Input  Processing. 

Input  processing  is  executed  once  per  frame. 


Input  data  is  screened  to  eliminate  known  fixed  targets  and  the 
baseline  returns  (if  desired).  This  is  accomplished  as  in  the  current 
system  by  assigning  gating  times  which  are  used  to  suppress  all  de- 
lay reports  in  the  desired  time  intervals. 

The  data  is  next  sorted  into  a set  of  input  tables,  one  for  each 
of  the  36  receiver/trarsmitter  combinations.  Each  data  word  con- 
sists of  a delay  time  and  a separate  bit  (normally  zero)  which  can  be 
set  to  indicate  that  the  data  has  been  utilized  for  track  processing. 


Figure  3.  Track  Files 

Track  information  is  kept  in  a set  of  track  files,  one  for  each 
receiver/transmitter  combination.  Each  file  has  room  for  a maxi- 
mum of  several  tracks  (four  tracks  have  been  alloted  in  the  initial 
program  version).  Each  track  slot  consists  of  three  words. 

The  first  track  word  (0)  contains  information  on  track  status. 

Bit  0 is  an  occupancy  indicator  showing  whether  or  not  the  track  slot 
is  being  used.  If  the  slot  is  not  occupied,  the  additional  information 
in  it  is  old  and  of  no  consequence.  Bits  1-2  indicate  the  status  type 
of  the  track  (see  Figure  8).  Bits  3-5  are  used  if  this  is  a parent 
track  (for  a track  split).  They  indicate  which  of  the  tracks  in  the  file 


is  the  trial  (split)  track  associated  with  this  parent.  Bits  6-10  hold 
the  firmness  of  the  track  (see  Figure  5).  Bits  11-15  contain  the  hit/ 
miss  count  which  is  utilized  in  track  initiation  and  drop  decisions 
(see  Figure  8). 

The  second  track  word  (1)  contains  the  smoothed  "position"  (actu- 
ally time  delay)  estimate  of  the  track.  It  is  the  quantity  used  directly 
in  vortex  position  calculations. 

The  third  track  word  (2)  contains  the  smoothed  velocity  estimate 
which  is  necessary  for  the  internal  operation  of  the  smoothing  algo- 
rithm. 

These  estimates  are  fixed  point  quantities,  fractions  whose 
scaling  can  be  fixed  once  and  for  all  by  considering  the  maximum  delay 
as  the  unit  of  position.  The  velocity  can  be  held  as  the  per  frame  change 
of  delay. 


Figure  4.  Smoothing  Equations 

The  smoothing  equations  recursively  update  position  and  velocity 
estimates  of  each  track  using  the  track's  data  for  each  frame.  The 
equations  are  a practical  implementation  of  a two-state  Kalman-Bucy 
filter. 

The  first  step  is  to  extrapolate  the  previous  smoothed  position, 
velocity  to  the  present  frame  time  (a  one-frame  advance).  Then  the 
deviation  between  the  datum  delay  time  dR  and  the  predicted  position 
(delay  time)  Dn  is  calculated,  and  fractions  (or,  0)  of  it  are  added  to 
the  predictions  to  produce  the  updated,  new  smoothed  position,  velocity. 

In  the  case  where  no  data  is  available  on  a frame,  the  extrapolated 
values  are  accepted  as  the  new  smoothed  values.  In  such  cases,  when 
a new  datum  is  finally  found  on  a later  frame,  the  8 constant  is  modi- 
fied by  dividing  by  the  number  of  frames  since  the  last  data  was 
received  (information  available  from  the  miss  counter).  Normally, 


The  selection  of  or,  0 smoothing  constants  is  implemented  by  a table 

lookup  which  is  keyed  upon  the  track  firmness  F . Firmness  ranges 

from  0 to  15  in  steps  of  1.  When  the  first  datum  is  received  for  a track, 

a firmness  of  0 is  assigned.  Subsequent  hits  increment  Fn  by  one, 

while  misses  decrement  it  by  two.  The  subsequent  F excursions  are 

n 

limited  to  the  range  1-15. 

As  the  firmness  increases,  the  appropriate  or^,  0 decrease  to  pro- 
vide increased  smoothing.  The  maximum  value  of  Fn  which  is  permitted 
is  a compromise  between  the  maximum  degree  of  smoothing  desired  vs. 
the  need  to  follow  quick  variations  in  the  delay  dynamics.  The  value  15 
is  only  an  initial  guess  of  the  appropriate  parameter. 

Search  bins  sizes  are  also  made  functions  of  firmness.  The  values 
listed  are  only  illustrative  and  may  be  calculated  from  the  known  stan- 
dard deviation  (a)  of  delay  noise.  The  bin  size  for  Fn  = 1 is  additionally 
based  on  the  expected  maximum  change  of  delay  position  per  frame.  Bin 
sizes,  &n>  may  require  experimental  adjustment,  since  the  theory  on 
which  they  are  based  oversimplifies  the  actual  track  dynamics. 


Figure  6.  a,  0 Values  Used  in  Delay  Tracker 


Figure  7.  Equivalent  Data  Weights 

The  a,  B smoothing  is  linear  for  constant  a,  0.  Thus  the  impulse 
response  of  the  filter,  the  weighting  of  the  past  data  values  to  achieve 
a position  (or  velocity)  estimate,  can  be  exhibited.  This  figure  shows 
the  position  estimate  weights  for  the  a,  0 at  a steady  firmness  of  15. 

It  reveals  that  the  weights  are  approximately  exponentially  tapered 
over  about  10  samples. 


Figure  8.  Track  Status  Transitions 

The  possible  track  statuses  are  tentative,  normal,  parent  and  trial. 

A tentative  track  is  created  when  data  is  initially  received.  It  re- 
quires N consecutive  hits  to  be  verified  as  a real  track  (rather  than 
noise).  The  status  then  changes  to  normal.  If  any  miss  occurs  during 
the  tentative  period,  the  track  is  dropped  and  must  be  reacquired  if  new 
data  arrives. 

A normal  track  is  maintained  until  too  many  consecutive  misses  are 
counted.  It  is  then  dropped  and  must  be  reacquired  if  data  continues. 

Under  special  circumstances  (no  data  in  primary  search  bin,  one 
datum  in  larger  secondary  bin),  the  normal  track  splits  into:  a parent 
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track  which  extrapolates  without  data,  and  a trial  track  which  follows 
the  data.  A special  or,  6 selection  may  be  appropriate  for  this  trial 
track  smoothing,  since  a fast  dynamic  variation  is  indicated,  and  of,  B 
should  not  be  too  small  to  follow  it. 

On  the  next  frame,  the  trial  track  survives  (becomes  normal,  and 
the  parent  is  dropped)  only  if  the  parent  receives  no  data  in  its  pri- 
mary bin  and  the  trial  receives  one  datum  in  its  primary  bin.  Other- 
wise, the  trial  track  is  dropped  and  the  parent  reverts  again  to  normal. 

Neither  parent,  trial  or  tentative,  tracks  can  be  split  in  this 
manner,  and  secondary  search  is  unnecessary  for  them. 

Figure  9.  Primary/Secondary  Correlation 

As  noted  above,  primary  search  is  applied  to  all  tracks,  and  the 
additional  secondary  search  is  used  only  with  normal  tracks. 

Bins  are  constructed  by  centering  at  the  predicted  position,  Dn> 
for  this  frame.  The  primary  tolerance  is  ± 6^,  and  the  secondary,  if 
needed,  is  ± 1.5  An-  The  enlargement  factor  can  be  experimentally 
adjusted. 
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Figures  10-11-12-13.  Correlation/Smoothing  for  Delay  Tracks 

These  figures  exhibit  the  proposed  cor  relation/smoothing  program 
flow  logic.  This  logic  is  entered  three  times  for  each  frame  process- 
ing of  a track  file.  Pass  1 processes  the  normal  and  parent  tracks  in 
a file.  Pass  2 processes  the  trial  tracks,  and  pass  3 processes  the 
tentative  tracks.  This  order  is  necessary  because  the  various  track 
categories  have  different  priorities  of  access  to  the  available  data. 

Data  which  are  used  in  any  pass  are  so  marked  in  the  input  file,  and 
cannot  be  used  by  a subsequent  track  or  in  a subsequent  pass. 

The  logic  contains  the  features  previously  described.  Some 
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additional  points  are  as  follows: 


a.  Before  changing  a tentative  track  to  normal, 
its  final  estimated  velocity  is  checked  against 
prescribed  limits  to  reject  all  tracks  with  im- 
possible velocities.  (Bin  size  also  limits 
excessive  velocity  excursions  implied  by 
the  first  two  data  points.  ) 


b.  If  a normal  track  has  more  than  one  datum  in 
the  primary  bin,  that  nearest  to  the  prediction 
is  selected.  However,  such  a condition  is 
counted  as  a partial  miss  (miss  increment  of 
1)  to  prevent  continued  tracking  through  heavy 
multiple  noise  reports.  (Ordinary  misses 
increment  the  miss  counter  by  2.) 

c.  If  a track  file  is  full,  splits  cannot  be  accommo- 
dated and  a normal  track  is  continued  as  a missed 
data  case  (not  shown). 

d.  Following  the  three  passes  through  primary / 
secondary  correlation,  some  unused  data  may 
remain  ,n  the  input  file.  This  data  is  utilized 
to  set  up  new  tentative  tracks  (initiation  pro- 
cessing). The  data  is  selected  in  order  of  the 
largest  delay  time  (as  is  done  in  the  current 
AVCO  system)  and  set  in  the  position  estimate 
word  of  an  unoccupied  track  slot.  The  velocity 
is  set  to  zero.  The  occupied  bit  and  other 
status  information  are  initialized  (status  = 
tentative,  1=1,  Hit  count  = 0).  Note  that  the 
insertion  of  position  data  is  equivalent  to  a 
smoothing  step  in  which  a=  1,  0 = 0. 


Background 

Figures  14-21  provide  some  background  material  relating  to  the 
firmness  control  of  a,  0 and  bin  sizes. 


Figure  14.  Variance  Factor  of  Position  Prediction  Errors  Due  to  Data 
Noise. 

This  figure  gives  contours  of  equal  position  prediction  smoothing 
in  the  a-0  plane.  Of  special  note  is  the  stable  triangle  of  or,  B values. 
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of  samples).  This  analysis  reveals  that  the  desirable  damping  factor 
of  Q = .7  produces  an  a,  g curve  close  to  the  previously  determined 
standard. 


Figure  17.  Optimal  a.  g Start-Up  Sequence 

The  foregoing  figures  and  analyses  were  based  on  steady-state 
assumptions.  By  utilizing  the  same  models,  the  problems  of  start-up 
and  missing  data  can  also  be  analyzed.  This  figure  shows  the  optimal 
start-up  sequence  of  a,  g ( small  circles)  which  begins  at  a = 1 , g = 1 and 
descends  asymptotically  to  a particular  steady-state  point.  By  repeat- 
ing this  calculation  with  various  relative  weights  on  noise/maneuver 
error  effects,  other  sequences  can  be  generated.  These  a,  B points 
all  lie  between  the  standard  curve  and  the  "start-up  envelope"  curve. 
Successive  points  fall  somewhere  on  the  heavy  traces  shown  in  the 
figure. 

Since  the  start-up  envelope  is  perturbed  only  very  modestly  from 
the  standard  curve,  we  adopt  values  of  a,  g on  the  latter  also  for  start- 
up. The  first  point  (Fn  = 1)  is  an  exception  where  the  true  start-up 
value  a=  1,  g = 1 is  utilized. 


Figures  18-19 

These  figures  illustrate  a similar  calculation  for  missed  data  in 
which  the  frames-since-last-data  factor  (Tn)  used  in  velocity  smooth- 
ing (see  Figure  4)  has  been  utilized.  Starting  from  the  steady-state 
condition  (SS),  it  successively  assumes  that  one  miss,  two  misses,  etc. 
occur.  During  the  misses,  the  track  is,  of  course,  extrapolated.  On 
the  first  datum  after  the  misses,  the  optimal  of,  g are  plotted  (labeled 
1M  after  one  miss,  2M  after  two  misses,  etc.  ).  On  the  second  datum, 
the  optimal  a,  g drop  to  near  the  steady-state  value. 

The  calculations  are  made  for  two  different  steady-state  starting 
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conditions  in  the  respective  figures. 

These  results  indicate  that  near  optimal  or,  0 for  missed  data  con- 
ditions can  be  selected  from  the  standard  curve  with  a suitable  firm- 
ness control  and  a T factor  correction. 
n 

Figures  20-21 

The  optimal  adjustment  of  bin  sizes  is  based  on  a comparison  of 
false  report  density  with  the  probability  density  of  the  predicted  posi- 
tion of  the  report.  The  width  of  the  latter  distribution  is  produced  by 
the  noise/maneuver  errors  described  previously.  If  data  is  missed, 
the  width  increases  but  the  height  of  the  distribution  decreases.  At 
first  this  increases  the  appropriate  bin,  but  finally,  the  bin  narrows 
and  vanishes. 

These  effects  are  built  into  practical  design  only  very  roughly. 

Thus,  the  enlarged  secondary  bin  is  utilized  only  in  a very  restricted 
way.  Also  the  vanishing  of  a bin  is  equivalent  to  the  decision  to  drop 
a track,  which  is  implemented  directly  by  counting  consecutive  misses. 

The  remaining  figure  illustrates  results  of  a calculation  of  devia- 
tion  (difference  between  data  and  predicted  position)  variance  during 
track  start-up.  The  different  curves  are  for  different  noise /maneuver 
weightings  and  lead  to  different  steady-state  values.  Since  the  suggested 
maximum  firmness  is  15,  the  appropriate  steady-state  parameter  is 
X».005  (see  Figure  15).  By  following  this  curve  and  multiplying  the 
standard  deviation  by  a safety  factor  of  3.5  to  4,  bin  sizes  listed  in 
Figure  5 are  obtained.  (Cf.  the  simpler  formula  given  in  the  descrip- 
tion of  this  figure.  This  formula  assumes  a X = 0 curve,  which  leads  to 
nearly  the  same  result.) 
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FIGURE  4. 
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FIGURE  5. 

TRACKING  PARAMETER  LOOKUP 


■ 

Firmness 

Position  Smooth 

Velocity  Smooth 

Search  Bin 

F„ 

an 

*n 

A 

n 

0 

1.000 

.000 

— 

1 

1.000 

1.000 

7.3 a v At  • v 

max 

2 

.833 

.700 

9.8a 

3 

.700 

.409 

7.3a 

4 

.600 

.270 

6.3a 

• 

• 

* 

• 

: 

• 

• 

• 

15 

.228 

.030 

4.6  a 

IF  HIT  n: 

Fn+l  = Fn+1 

(F 

n max 

IF  MISS: 

Fn+l=Fn'2 

(F 

' n min 

INITIAL  HIT: 

F = 0 -*  F,  = 1 

O i 

B-19 


16 


EQUIVALENT  DATA  WEIGHTS  FOR  SMOOTHING  AT  MAXIMUM  FIRMNESS 
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TRACK  STATUS  TRANSITIONS 
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FIGURE  9. 

PRIMARY/SECONDARY  CORRELATION 

• PRIMARY  SEARCH  - ALL  TRACKS 

• SECONDARY  SEARCH  - NORMAL  TRACKS 
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FIGURE  17. 

OPTIMAL  a,  /?  START-UP  SEQUENCE 
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FIGURE  20. 
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APPENDIX  A 


Delay  Tracking  Program 

This  program  was  designed  to  improve  system  performance  by 
enhancing  the  reliability  and  accuracy  of  delay  tracking.  The  im- 
provement results  in  an  overall  increased  immunity  to  noise  and  a 


reduced  susceptibility  to  data  discontinuities.  This  program  is  a 
simplified  version  of  the  complete  system  described  in  the  body  of 
the  report. 


The  program  consists  of  six  routines: 


1.  Main  routine  ARCON  provides  the  bulk  of 
logical  and  control  operations  of  all  other 
routines. 

2.  Routine  INITIA  performs  the  process  of  new 
track  initiation. 

3.  Routine  INSERT  inserts  new  tracks  which  are 
initiated  by  INITIA -routine  into  track  file. 

4.  Routines  LOOKTR  and  QUEST  look  for  a track 
of  the  proper  status  which  should  be  con- 
tinued and  the  data  for  the  job. 

5.  Routine  DROP  eliminates  tracks  which  cannot 
be  continued  properly  (see  flow  chart  of  the 
algorithm). 
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OS/360  FORTRAN 


COMPILER  OPTICNS  - NAM  E = HA  I N, OPT  = 0 1 , L I NECNT*50 , S I ZF=OOCDK , 

Snua.CE ,.EBCD IC., NOLIST , MOQE.CK. .LOAD  , MAJS, N3E0 mjfl, NO XRE£_ 

> SUBROUTINE  ARCON( KDAT  A, NT  RACK ,NMB»NFRM, NMBR .KOOTR.P ARAM.9T ) 

I_ DIMENSION  TRC.1 18,  A»2),NST(18»A,5)  . 

V DIMENSION  DTI18,ll),NTRACKI2)»K3DTR(2),NMRO  (18, A,2> ,KOAT4 (18) 

i ; OIMEKS.IOM..P.ARA.H.I L5.3)  ....  _ 

i NCATA=10 

K T 


* NHITLM=*3 

i VJ£AX=1C. ...J.;... __ 

> 1F(NF«M)1, 2, 1 

L 2 MER.M-N  FRME.1 1 ! - 

\ OC  7900  NN- 1*18 

I EC  fc 700  KK » 1 . 4 

i DC  6600  LM=1,5 

>  66.Q0UMSTJMN « KK.*  LM).*0 

i 6700  CONTINUE 

i iflao-ccwiinuE . 

1 CALL  INITIA<OT,TRC,NST,NTRACK,NOATA,NT,NMBR ,NMB,NFRM) 

2  QG-.ll.JLg-1-fcIB 

? KCATAin-1 

I LI Bill.  LL-77.77. . 

J RETURN 

I 1 H£BM«NEP,M«-.l 

t CC  3 NPASS 3 1 »2 

5 L0EG.LN.-_1 '. 

f;  J»ECINxO 

I l.-NTRACKI  NPASS ) ■_ 

i JMIU.3.4 

1 » NPT-KODTRl  NPASS  1 . ..  . _ 

J CO  5 KP-1,1  ’ ' 

I Rfi.U.E.(6, 5QQ.J  N.TR  ACX 

1 500  FCRMATI2I10) 

L HAm.Jl  6,3.5  IlBEGINl.  JBEGIN,  I.NP.T. 

t,  55  FCRMATIIX.AIIO) 

5 CAL  L_L.CGKI_R.LI  BEG.  I N » J BFG I N_«JNP  T. , N S T , N T.L  

S P«ED*TRC( I BEGIN, JBEGIN, II+TRCI I BEGIN, JBECIN.2  » 

Z U£ILUSiP£EQ±P-A^AtJ  INST  J.  I B_E£  JN , J BE.G.I.N*  .3.) .,  3J 

[)  OCWBN-PREO-PARAM(NST(  IBEG  IN,  JBEGIN,  3)  , 3 I 

2  CALI QUEST  .1.1. BEG  IA  , DT,  MOAT  A , DOM  BN  ,UPBIN,  PRFD  , I.Q,  N.OI 

) 1F(N0-1)6,7,B 

L 6 1 F (NPASS-1  i 3 ill,  1 3. 

? 5 IF(NP ASS-113, 9, 13 

3  9 \ST(  1 BEGIN  .JBEGIN,  A ) *NST(  I BEG  IN . J BEGI N . A ) *1 

4 I FINS  TUBE  GIN,  JBEGIN,  A ) -MSL  IM)18, 18,13 

5  L2. N S.TJ.I.B.EG LN_.  JBEGIN  , A ) *NST. I TBEGIN,  JBEGIN,  A )+2 

S IF( NS T(I BEGIN, JBEGIN, A) -MSL IMIlOtlG, 13  - 

7 IQ LE.INSLI  I.BEGIN,  J..BEGIN,  3J-2  ) IA,  1A,  15.  ■ 


I 


48 

.49 

14 

NSTII3EGIN, JBEGIN, 3)=1 

_G.C_LO  5 _ 

5C 

15 

NSTI I BEGIN, JBEGIN, 3) =NST( I BEGIN, JBEGIN,  3 )-2 

—GJtL.T.G  5 _ 

52 

13 

NST ( IBEGIN, JBEGIN, l)*0 

53  NT.RACMNP.ASS  ).«NTRACX<  NPA$S)-1 ...  _ _ 

54  CALL  ORCP(NST ,NMBR, IBEGIN, JBEGIN I 

.55 fi.C-IQ._5 

56  7 IF(NPA$S-1)3,16,17 

57  15 fi.SK  IBEGIN,  JBEGIN.,  4)«0  „ 

58  18  TRCIIBEGIN, JBEGIN, 1)=PRED4PARAMINST(IBEGIN, JBEGIN, 31,1)* 

JUUtT.<m6IN..I01-PRfO|  

55  TN=NST(  IBEGIN,  JBEGIN,  4H-1 

6C L8LLiaf-Gl.N-.jnEG IN , 2 )=TRC(  rBEGINf.JJ3fcG.tH,.2  ft. 

1P«RAM(NST( IBEGIN, JBEGIN, 3) ,2 )*(DTC l BEGIN, IP)-PRED)/TN 

61  1F.1NP  ASS.-1.13 ,24,19.  ..  ... _ _ 

62  17  NST( IBEGIN, JBEGIN, 5)=NST( IBEGIN, JBEGIN, 5 M-l 

63  GC._IQ  .ia. ____ 

64  19  IFINSTUBEGIN,  JBEGIN,  5 1-NH 1TLM  1 24, 27, 27 

65. 2.7. I EJ  TR  C.UR.EG.I  N...JBE  G IN, .2  )1 20 , 21,71 

66  20  V= ( -TRC ( IBEGIN, JBEGIN, 2) ) 

6.7. JLC.TQ.  22..* ..  . ...  . ..  ...  

68  21  V=TRC< IBEGIN, JBEGIN, 21 

6.5.  22 .LE.IV-VMAXI  23,23, 13  _ 

70  23  N ST (IBEGIN , JBEG IN ,21*2 

71  N5KI BECLLN , JftEGIN,  4)  = C! 

72  N TRACK  C 1 ) = NTRACM  1)41 

73. IL1RACK(2)  = NTRACK(  2 )-l  ...  _ 

74  24  IMNSK  IBEGIN, JBEGIN, 31-15)25,26, 2S 

7.5.  75 N.SIUB.EGJ.N, . JBEGIN,  3)=NST(  IBEGIN,  JBEGIN,  3)41  _ 

76  CALL  INSERKTR.C, NST, NM9R,NFRM, IBEGIN, JBEGIN) 

.27 25  PC  2B  KA*  I P. NDAT  A 

78  0T(IB6GIN,NA)*t)T(  IBEGIN,NA*1I 

79  I F (.07  I IBEGIN, NA  1-7777,  >28,  5,28  . _ 

so  23  continue 

.81...  5 CCNTI.NUE 

B2  3 CONTINUE 

83 C AJ.L ..  JN ITJA I UT  , TRC , NST,  NTRACK , N D A T A » N.  T,  N MBR  » NHB,  NjFJjt  M ) 

e4  CG  30  1 = 1,  18 

85  .KCM.A(I)  = 1 _.A.  _ 

86  30  C T( I • 1 1*7777. 

F7_ RETURN  _ . _ . ..... 

88  END 
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— tFIDTI I .L1-77Z7. ) 13 

28 

13 

CCNTINUE 

20 

8 - 

CfLNtl.NUF 

30 

18 

I F CMH ) 1 , 15  * 19 

XI 

i q 
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16 
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15 
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36 
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37 
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38 

A ST  1 1, L, 5 1=0 
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AO 

TRC(I,L,1)=DIN(J) 

A1 
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OS/36)  FORTRAN  H 


COMPILER  OPTIONS  - NAME*  MAIN,0PT*01 tL INECNT-50 ,S l ZE-OOOOK , 
SUllRCF  >E  ftf  Off  «NOL  1 ST  ►.NTOfiCK.  »l  flAO  P » NTIEQ  LJ  » IOtNflXRF_F 


C2 

C3— 

SUBROUTINE  LOOKTRI IBEGIN, JBEGIN.NPT.NST ,NT) 
riMF^MON  N9T  I \Rt  *.,«>> 

04 

C5 

5 

IF(JBEGIN-NT)5,'a,1 
/ncr.jMj.  j^pr.  |n+i 

C6 

CJ 

4 

GC  TO  7 

IRFr.lKsI 

C 8 

I SEC  IN* IBEG!N*1 

(15  - 

...  X- 

nr  j ts_iBFGiNt.ia 

1C 

CC  2 J* JBEGINi NT 

LI 

I F INST  1 1 t J» 1 1 1 1 > ? > 3 . . _ 

12 

3 

IP(NST(I*j*2)“NPTj2fA»2 

1 ■* 

? 

1A 

JBEGIN*1 

L5l 

r rjNT  INMP 

16 

A 

IPEGIN*I 

11 

pppr.iN».i 

18 

RETURN 

1.9 — 

_EJiO.  
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CONTINUE 


COMPILER  OPTIONS  - NAME*  MAIN.OPT-Ol tL INECNT-50 . S ! 2E»OOOOK , 

SOURCE.  EBCDIC. NO!  I ST . NnnPf.lt  . I DAD  ■ MA P . NOED I T . t n. NOXRCF 

'02  SUBROUTINE  INSERT  I TRC • NST , NMUR, NFRM, I , J ) 

03 DIMENSION  TRC  MB  .A.21.NSTI  1B._4.5-1  .MHBRIlB.A^l) 

OA  C CMNON/AAA/ IFRM(lB,4,2G)t RESULT I IB. 4, 20) 

C5 ‘'FIRM*  5 

06  NMBRt  l,J,2)*NMBR(  I.J.2H-1 

C2 1 FRHC  1 . J .MMBRI 1 « J « 211  »NF RM 

CB  RESULT (I .J.NNaRI I, J, 2)|>TRC( I.J.ll 

10  ~ 2 !FINSTII,J,3)-MFIRM)3,4,4 

11  4 WRITE  16.10  I NM3 RJ.Lt  J_»l_l  ,_LI £RH  l ±%J-ULLt HsljJLQ,  JL) 

12  10  FORMAT 1 1 10, 101 10) 

13  HRITFI4.il  I RESULT  lI.J.ll.lRESULIll.J.ILl,  K*2.»  10,  1 1 

l A 11  FGRMATIF20.2.9F10.2) 

1C  U R f TP  14. 1711 FRM (I.. 1.111. II ERM ( I.J.KI.K: 12.20. 1) 
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COMPILER 

OPTIONS  - NAME"  HA l N» 0PT«01 »L 1 NECNT  = 50»  S I ZE*0003K  » 

‘ ' 

)C2 

-C3 

SUBROUTINE  CROP (MS T, NMBR , I , J » 

OIW£NS2  QN_.NSI.L18 4 ,5J.,NMBRlia  ,4.71 

H'4 

105 

CCMMON/AAAZIFRM(18,4,20),RE$ULT<ia,4,20) 

-KF-IRMg.5 

>C6 

<2 

_J 

I FtNSTt I »J,3)-MFIRM)2,1»1 

IFINMRRII.J. 21-1013.3. 4 

■ C8 

C9 : 

3 

N=KM8R  ( I , J , 2 ) 

HJ1LI.E..(j6  «.1Q  J.NMBJLL.L,.  JL»  LL,_l  LF.RM  11j_J  . K)  ,K=L.M.l) 

10 

1J. 

10 

FCRHATU  10,10110) 

— HEJXEX6xU-lE.ESULTJX».Jf  1 It 1AESUUXL»JjKJ_.J1h2  .H,JLL_ 

12 

13 

11 

FORMAT ( F20 .2,9F10.2) 

GO  TO  7 

14 

15  

4 

WRITE (6,10 )NMBR( I ,J,1 ),( IFRMt I,J,K) ,K-1 ,10,1) 
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1J 

WRITE (6, 11 JRESUL  T ( I , J , I ) , IRESULTC I, J,K) ,K=2,10,1) 

— W JLLI6.(A,.12LLLERN  l LtJ.*  LL)  JL2jM jJLl 

18 

13 

12 

FORMAT ( 120,9110) 

WRITE (6,11 IRESULT1 1,J,11), (RESULT 11 ,J .K ) .K-12.M.1) 

20 

2.1 

2 

N MRR( I ,J,2)»0 
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22 
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